Orbitally induced string formation in the spin-orbital polarons 
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We study the spectral function of a single hole doped into the ab plane of the Mott insulator 
LaVOa, with antiferromagnetic (AF) spin order of S = 1 spins accompanied by alternating orbital 
(AO) order of active {dy,,dz^} orbitals. Starting from the respective t-J model, with spin-orbital 
superexchange and effective three-site hopping terms, we derive the polaron Hamiltonian and show 
that a hole couples simultaneously to the collective excitations of the AF/AO phase, magnons and 
orbitons. Next, we solve this polaron problem using the self-consistent Born approximation and 
find a stable quasiparticle solution — a spin-orbital polaron. We show that the spin-orbital polaron 
resembles the orbital polaron found in eg systems, as e.g. in K2CUF4 or (to some extent) in LaMnOs, 
and that the hole may be seen as confined in a string-like potential. However, the spins also play a 
crucial role in the formation of this polaron — we explain how the orbital degrees of freedom: (i) 
confine the spin dynamics acting on the hole as the classical Ising spins, and (ii) generate the string 
potential which is of the joint spin-orbital character. Finally, we discuss the impact of the results 
presented here on the understanding of the phase diagrams of the lightly doped cubic vanadates. 

PACS numbers: 71.10.Fd, 72.10.Di, 72.80.Ga, 79.60.-i 



I. INTRODUCTION 

A single hole doped into the half-filled Mott insulat- 
ing ground state with antiferromagnetic (AF) order does 
not move freely as its motion disturbs the spin order ii In- 
stead, it couples to the collective (delocalized) excitations 
of the AF ordered phase (magnons), and it propagates 
through the lattice together with a "cloud" of magnonsi^ 
Thereby the energy scale of the "coherent" hole propa- 
gation is strongly renormalized and is given by the AF 
superexchange constant J. In this way a quasiparticle is 
formed which is frequently called in the literature a spin 
polaron^ Actually, such a phenomenon is not only of the- 
oretical importance but also describes a realistic situation 
found in the photoemission spectra of the parent com- 
pounds of the high- Tc cuprates, as e.g. in Sr2Cu02Cl2i^ 

A more complex situation can occur in the systems 
with partly filled degenerate orbitals, where a doped hole 
may not only couple to magnons but may also couple to 
the crystal-field excitations 1^ Then we may for example 
arrive at a problem shown schematically in Fig. [l] which 
is studied in this paper and discussed in more detail be- 
low, or end up in two somewhat simpler, though still 
challenging, situations where at least the coupling to the 
magnons could be neglected. To get a better insight into 
the orbital physics we give a brief overview of these two 
latter cases in the next two paragraphs. 

First, we recall an example of the ab planes of LaMnOs 
which have ferromagnetic (FM) spin order and alternat- 
ing orbital (AO) order of Cg orbitals in the ground statci^ 
It has been shownl that although the hole introduced 
into such a state does not disturb the FM spin order, 
it couples to the collective excitations of the AO state 
(orbitons). Here again a quasiparticle is formed which is 
called this time an orbital polaron^ It should be noted, 
however, that the orbital polaron has an even smaller 



bandwidth than the spin polaron^i as the orbitons are 
in general much less mobile than the magnons (or even 
immobile) due to the lack of the SU(2) symmetry in 
the orbital systems;^ and almost directional Ising-like 
superexchangci^ii^ Actually, one can understand the hole 
motion in this case in terms of the string picturei^i The 
hardly mobile orbitons cannot repair the string of the 
misaligned orbitals in the AO state, which arises by the 
hole motion. Thus, it is the hole which has to return 
to the original site and cure the defects by retracing its 
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FIG. 1: (Color online) Artist's view of a single hole introduced 
into the spin and orbitally ordered ab plane of LaVOa which 
leads to the spin-orbital polaron considered in the paper. The 
electrons occupy the dyz and d-x degenerate orbitals forming 
the classical AO state (their projections along the a and 6 
axis are shown) whereas the electron spins alternate on the 
neighboring sites forming the classical AF state. Only the spin 
1/2 of the electrons in the dyz and dzx orbitals are shown. At 
each site an additional electron occupies the d^y orbital (not 
shown). Due to Hund's coupling the two electron spins form 
a total spin 5 = 1, while at the site of the hole S = 1/2. 
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path, unless it propagated due to small off-diagonal hop- 
ping in an Cg system^ and no defects were created on 
its path (the latter process also contributes to the above 
mentioned very small bandwidth of the orbital polaron) . 

Second, only recently an even more extreme situation 
of the system with orbital order was investigated i^^i^^ 
The case of a hole doped into the plane with FM spin 
order accompanied by the t2g AO order (which could cor- 
respond not only to the hole introduced into the ordered 
ground state of Sr2V04 with t2g orbitals but also, sur- 
prisingly, to those of K2CUF4 or Cs2Ag4 with Cg active 
orbitals, see Ref. [Tsh . Also here a quasiparticle is formed 
due to the dressing of a hole by the collective excitations 
of the ground state with AO order. However, due to the 
specific t2g orbital symmetries the orbitons are not mo- 
bile at all, the off-diagonal hopping is prohibited, and 
the quasiparticle acquires a finite bandwidth only due 
to the frequently neglected three-site terms.— Thus, the 
string picture dominates the character of the t2g orbital 
polarons even more than in the case of systems with eg 
orbital degrees of freedom. 

Finally, after this brief overview, we turn to the prob- 
lem investigated in the present paper: the motion of a 
single hole introduced into the ground state with coex- 
isting AF order of = 1 spins and t2g AO order, shown 
schematically in Fig. [1] This situation corresponds to a 
single hole doped into the ab plane of LaVOaiii It dif- 
fers from all of the three (one — purely spin and the 
other two — purely orbital) cases described above due 
to the coexisting AF/AO order. Neither spin nor the or- 
bital background is then transparent for the propagating 
hole, and the hole has to couple both to the magnons and 
orbitons simultaneously (see Sec. Ill Cp . To our knowl- 
edge, the mechanism and physical consequences of this 
coupling have never been studied beforeii^ 

Furthermore, the coexistence of the AO and AF or- 
der is an extremely rare case in nature as it formally 
violate&i^ the Goodenough-Kanamori rules that state 
complementary spin and orbital order in the ground 
state, i.e., either FM spin coexisting with AO order, or 
AF spin coexisting with FO order — thus, it is interest- 
ing to investigate how a hole can propagate in a state 
which does not fulfil this rule and by its nature contains 
more quantum fluctuations. Therefore, the main goals of 
this paper arc to verify: (i) under which conditions, if at 
all, a quasiparticle can be formed in such an AF phase 
with AO order, (ii) if a quasiparticle is formed, w/iai are 
its properties and whether the string picture (which is 
typical for orbital physics) can explain them, and finally 
(iii) what is the role of spins and orbitals in the possi- 
ble formation of the string in the spin-orbital polarons. 
By working out the answers to the above questions we 
would like to predict the main features of the photoe- 
mission spectra of the cleaved samples of LaVOa with 
polarization corresponding to the ab planes. 

Before we move on to present the answers to these 
questions, let us note another motivation for studying 
the above problem. It concerns the phase diagram of 



the lightly doped cubic vanadate Lai-xSrxVOsi^ It ap- 
pears that in this strongly correlated compound the C- 
AF (AF ab planes with FM order along the c direction) 
coexisting with G-AO (AO order in all three cubic direc- 
tions) order Mott insulating phase is not only stable at 
X — but also persists to a relatively high value of dop- 
ing X = 0.1781^^ Furthermore, the C-AF phase remains 
stable up to an even higher value of a; = 0.26, although in 
this regime the insulating and orbital ordered phase has 
already disappearediii As in the ionic picture x stands 
for the fraction of holes doped into the d orbitals of vana- 
dium ions (where a nominal valence upon doping changes 
as d^~^), it remains a challenge to explain why the or- 
dered and insulating states persist to such high doping 
level Besides, somewhat similar phase diagrams have 
been observed in other doped cubic vanadates such 
as Pri_^Caj;V03, Ndi_a;Srj;V03 or even to some extent 
in Yi-^i-Ca^i-VOa although in all these cases the lattice 
(Jahn- Teller and GdFeOa-like) distortions contribute sig- 
nificantly to their physical properties, similar to the un- 
doped parent RVO^ (with i?=Pr,Nd,Y) compoundsi^ 
Here to avoid further complications due to additional or- 
bital interactions imposed by the lattice, we would like 
to concentrate ourselves on the spin-orbital polarons in 
(almost) undistorted compound LaVOs. 

Actually, one should in principle be able to understand 
some of the features observed in the phase diagram of the 
doped cubic vanadates by comparing them with those 
of the high-Tc cuprates or the colossal magnetoresistive 
manganitcs. However, such a comparison only further 
emphasizes the lack of theoretical understanding of the 
doped cubic vanadates. First, in the high- Tc cuprates, 
such as La2_a;Sra;Cu04, the AF order disappears very 
quickly with doping x, i.e., already for x ~ 0.023i This 
is despite the fact that the value of the superexchange 
constant J is relatively high there which naturally leads 
to larger magnetic energy in the cuprates than in the 
vanadates. This suggests that, among other factors, it 
is the orbital dynamics which could be responsible for 
the totally different behavior of these two classes of com- 
pounds upon hole doping. Second, a similar conjecture 
can be drawn from the comparison between the vana- 
dates and the manganitcs. In the latter ones, e.g. in 
Lai_2;Sr2;Mn03, the AO orbital Mott insulating state is 
stable up to a; 0.18;^ i.e., almost to the same doping 
level as in Lai-^jSr^jVOa. However, Lai-a^Sr^^MnOs has 
FM planes but Lai_3;Sr2;V03 has AF planes - this again 
indicates that it is the orbital dynamics which governs 
the behavior of holes in the doped cubic vanadates. 

Hence, we arrive at the following paradox. On the one 
hand, it is obvious from the theory that any meaning- 
ful study of the doped pcrovskite vanadates should take 
into account the spin and orbital degrees of freedom on 
equal footing. On the other hand, we see from the above 
conjecture from the experiment that these are the or- 
bital degrees which are, to large extent, responsible for 
the stability of the observed phases of the lightly doped 
cubic vanadates. Thus, another aim of the paper is to 
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shed some light on this apparent paradox. Certainly, 
the studies presented below concern only the case of just 
a single hole doped into the half-filled compound and, 
furthermore, they investigate only the situation in the 
two-dimensional (2D) model focusing on the ab planes. 
Nevertheless, we hope to understand the generic role 
of the spin and orbital degrees of freedom in the hole 
doped cubic vanadates. Besides, we note that comple- 
mentary studies presented in Ref. revealed the role of 
the AO and FM order stable along the third (not stud- 
ied here) direction in the doped Lai-^jSra-VOs, and ex- 
plained the differences between the doped Yi_a;Caa;V03 
and Lai_xSr2;V03 but, by the very nature of that one- 
dimensional model, could not address the problems men- 
tioned above. We also stress that, contrary to the prob- 
lem solved in this paper, the hole moving along this 
third c cubic direction couples only to one type of exci- 
tations: either orbitons in the lightly doped C-AF phase 
of La1_2.Sr2.VO3, or magnons in the very lightly doped 
G-AF phase of Yi_^Ca:rV03,^ 

After a thorough discussion of the motivation standing 
behind the studies presented in the paper let us now con- 
centrate on the model and the theoretical method used in 
this work. We use the 2D t-J model for the ions in 
the configuration with degenerate t2g orbitals, and the 
superexchange given in Refs. [Hill, as the model which 
describes the physical properties of a hole in the ab plane 
of the lightly doped cubic vanadates. Since any t-J model 
which contains Ising-type interactions should be supple- 
mented by the three-site effective hopping terms in order 
to provide a possibility of coherent hole propagation^i^iH 
we derive these terms and add them to the t-J model 
(such a model is then called the strong-coupling model). 
Next, we solve the strong-coupling model in the case of 
a single hole doped into the ground state of a Mott insu- 
lator by reducing it to the polaron-like model and intro- 
ducing the self-consistent Born approximation (SCBA)i^ 
Although the final solutions (spectral functions) are ob- 
tained numerically on a finite mesh of the momentum k 
points, the SCBA method is known to be very accurate 
to study such problems and is largely independent 
of the size of the clusters^ on which the calculations are 
performed. 

The paper is organized as follows. We start with in- 
troducing the spin-orbital t-J model in Sec. Ill Al which 
is supplemented by the crucial three-site terms for the 
spin-orbital model, see Appendices [X] and [B] for detailed 
discussion and derivation. We reduce the derived strong- 
coupling model {t-J model with three-site terms) to the 
spin-orbital polaron problem in Sees. Ill Bl - InD] and dis- 
cuss its parameters in Sec. Ill El Next, we formulate the 
solution of the model using the analytic approach within 
the SCBA method, see Sec. IIIII In Sec. IIVI we present 
the numerical results and discuss the properties of the 
quasiparticle states. Finally, in Sec. |V]we analyze the 
composite interplay of spin and orbital degrees of free- 
dom in the formation of the spin-orbital polaron whereas 
in Sec. IVII we present our conclusions. We adopt a con- 



vention throughout the paper [despite the introductory 
part (Sec. [I]) and Appendix [K\ . that the red (blue) color 
in the figures denotes features associated solely with spins 
(orbitals). 



II. THE MODEL 

A. The spin-orbital t-J model with three-site terms 

The Hamiltonian of the spin-orbital t-J model with 
the three-site terms (called also the strong-coupling spin- 
orbital Hamiltonian), relevant for the planes of lightly 
doped cubic vanadates, consists of three terms, 

H = Ht + Hj + H3s, (2.1) 

where the last part H^s stands for the three-site terms. 
The first term in Eq. (|2.ip . which describes the hopping 
of the electrons in the constrained Hilbert space, i.e., in 
the space with singly occupied or doubly occupied sites, 
follows in a straightforward way from the unconstrained 
hopping of electrons residing in the t2g orbitals. The lat- 
ter one has to fulfill the t2g orbital symmetries meaning 
that along each bond in the cubic direction only two out 
of three t2g orbitals are activeJ^i^^ This means that in 
the ab plane, which is under consideration here, electrons 
in dyz = a {d^x = b) orbitals can hop only along the b (a) 
direction. Besides, the dxy orbital does not contribute 
to hopping elements as it lies lower in energy and is al- 
ways occupied by one electron in the half-filled and lightly 
hole-doped regime;^ Hence, we arrive at the kinetic oc t 
part of Hamiltonian (|2.ip , 

Ht = -tJ2'P {blbi+^.^ + «L«i+b,. + H.c.) r. (2.2) 

i,cr 

Here the constrained operators 

bl = C(l-"ifc5)(l-"iaff)(l""iaa), (2.3) 

«L = aL(l -"■;"*)(! -"ibff)(l-»ifca), (2.4) 

mean that the hopping is allowed only in the restricted 
Hilbert space with not more than one {a, 6} electron at 
each site i (a stands for the spin component opposite to 
a). Besides, since the Hund's coupling is large {Jh ^ 
t) in the cubic vanadates;^ we project the final states 
resulting from the electron hopping onto the high spin 
states, which is denoted by the V operators in Eq. (|2.2p . 

The middle term in Eq. (|2.ip describes the spin-orbital 
superexchange in cubic vanadatesi^ — here we use the 
notation of Ref. [H, see Eqs. (6.5)-(6.7). It reads 

H,j = h[P +Hf +Hf\ (2.5) 
where the contribution due to different excitations in the 
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multiplct structure of vanadium ions arei^ 



motion read: 



H 



(1 



-ijri5:(Si.Sj+2)(i-71^r/), (2.6) 

(ij> ^ ^ 



(ij> 



(2.7) 



H 



(3) 



Jr3 (Si • Sj - 1) (I T T + 2T^T, 
^ (ij> ^ 



(2.8) 



Here: Si is a spin 5 = 1 operator, Tf = (riit, — nia)/2 
is a pseudospin T = 1/2 operator, and J ~ it'^/U is the 
superexchange constant, with U being the intraorbital 
Coulomb element on a vanadium ion. As in the spin 
t-J model)2£ we assume that the hopping t is small, i.e., 
t <^U whose assumption is well satisfied by the realistic 
parameters in the cubic perovskitesiii The factors ri — 
1/(1 — 3?7) and = 1/(1 + 2ri), where 



V = 



U 



(2.9) 



account for the Hund's coupling Jh and originate 
from the multiplet structure of the excited states 
characterized by the various allowed spin and orbital 
configurations!^ Let us note, that superexchange Hamil- 
tonian (|2.5p was derived^! with the assumption that the 
dxy orbital was singly occupied at each Vanadium ion 
(see discussion above). Besides, in principle Hamiltonian 
(|2.5p was originally derived for the undoped case and 
should be modified for the doped case by adding the su- 
perexchange interactions due to the existence of the d\S^ 
and djdj nearest neighbor configurations. However, in 
the discussed here small doping regime, and particularly 
for a single hole, the modification of magnon and orbiton 
excitations due to these terms should be very small and 
we will neglect them. 

The last term in Eq. (|2.ip denotes the three-site terms 
which, to our knowledge, have not been derived before for 
the case of the systems with spin and orbital degrees 
of freedom. These terms, although frequently neglected, 
are crucial for a faithful representation of the Hubbard 
model in the strong-coupling regime;^ and can play an 
important role in the coherent hole propagation in orbital 
systems J^ii^ However, in the present case the derivation 
of all possible three-site terms is relatively tedious and 
leads to a complex expression. Fortunately, it was shown 
in Ref. [l^ (for orbital systems) and in Ref. (for spin 
systems) that the only three-site terms which occurred 
to be relevant for the lightly doped systems were these 
which did not contribute to the coupling between hole 
and orbital or spin excitation but merely contributed to 
the free hole motion (see Sec. IH D[) . Nevertheless, we 
discuss the possible role of the neglected three-site terms 
in Appendix |B1 Then, the derived in Appendix [X] three- 
site terms which could possibly contribute to the free hole 



\.<7 

(2.10) 



Note that these terms are oc J and hence are of the same 
order in t'^ /U as superexchange terms p.5|) . 

In the next four sections we reduce the model Eq. (j2.ip 
to the polaron problem following Ref. 0. This would 
enable us to solve the model Eq. (12. 1|) using the SCBA 
method. 



B. Undoped case: low energy excitations 

It was shown in Ref. [3 that the ground state of the 3D 
spin-orbital superexchange model (|2.6p is characterized 
by the classical phase with C-AF spin order and G-AO 
orbital order. Thus, we start from the above undoped 
ground state of the Hamiltonian Eq. (|2.1|) with classical 
(Neel) order. As usual, this is not the eigenstate of the 
Hamiltonian and thus the full description of the system 
has to take into account the quantum fiuctuations around 
such a classical ground state. Below, we will calculate 
them, together with the low energy (magnon and orbiton) 
excited states, employing the linear spin wave (LSW) and 
linear orbital wave (LOW) approximation, following Ref. 

First, in the classical state we introduce two sublattices 
A and B such that: (i) all a (b) orbitals are occupied in 
the perfect AO state in sublattice A (B), and (ii) spins 
pointing "downwards" ("upwards") are located on sub- 
lattice A (B), as shown schematically in Fig. [1] Next we 
rotate spins and orbitals (pseudospins) on sublattice A 
so that all the spins and pseudospins are aligned in the 
considered ab plane and take positive values now. 

Second, we introduce Schwinger bosons t and / for 
orbitals and spins such that: 



^iflh^-flfil): 



with the constraints 



E 

7— a,b 



1 , 



^ y /io-/i£ 



(2.11) 

(2.12) 

(2.13) 
(2.14) 

(2.15) 



Third, we transform the Schwinger boson operators into 
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the Holstcin-Primakoff bosons a and /3: 



4 = \/l-^Liia^\/l-A% (2-16) 



4 + = 



(2.17) 



/iT = V2-//i/u^ V2-"l"i' (2-18) 
fl=al (2.19) 

where the above constraints arc now no longer needed. 

Next, wc substitute the above transformations into the 
Hamihonian Hj and skip higher order terms (using LSW 
and LOW approximation). The latter approximation 
physically means that the number of bosons ai and 
which describe the fluctuations around the ordered state, 
is small. This results in the effective substitutions: 



1 



- 2 ~ 
S? = 1- alai , 

5;+ = V2ai , 

Sr ^ V2al . 



(2.20) 

(2.21) 
(2.22) 
(2.23) 



Finally, we introduce Fourier transformation sepa- 
rately for each sublattice {N is the total number of sites 
on both sublattices and the lattice constant a — 1): 



(2.24) 
(2.25) 
(2.26) 



j6B 



we define operators 

ak± = ("kA ± aks)/\/2, 



(2.28) 



and perform the standard Bogoliubov transformation to 
diagonalize the magnon part of the Hamiltonian: 



fiki = -"kaki - I'kalki 



(2.29) 



where 



Wk = \/ ^t^^ > = ~sgn(7k)W '^^ , (2.30) 



2i/k 



2i/k 



with i/k = -^/l — 7^ and the structure factor for the 
square lattice 7k = (cos + cos ky)/2. 



Then, after neglecting constant terms, the LSW and 
LOW Hamiltonian for magnons and orbitons reads: 

'Hj = Hi^sw + ^Low, (2-31) 

-f^Lsw = '^k(ak+Q!k+ + al-ttk- + 1)> (2-32) 

k 

i^.ow = Jo Y.^PIaP^.a + PIbP^b): (2.33) 

k 

where 

~ " (2.34) 

and 



Js 



' (l-377)(l + 277) 

1 - 377 - 577^ 
4(l- 377)(l + 2,7) 



J. 



(2.35) 



in agreement with Eq. (6.11) of Ref. l24| and Eq. (11) of 
Ref. [2^. Besides the dispersion relation of the magnons 



(2.36) 



where z = 4 is the coordination number and S" = 1 is the 
spin value. Let us note, that in the regime of reasonable 
values of ?7 g [0,0.20]: Jo > and Js > 0, which means 
that the classical ground state has indeed coexisting AF 
and AO order. Furthermore, at temperature T = the 
considered here classical 2D AO and AF ground state is 
stable with respect to the quantum fluctuations, both in 
spin and in orbital channel. In fact, the orbital order is 
undisturbed by local Ising excitations, while the quantum 
AF ground state is modified and the order parameter is 
renormalized with magnon excitationsi^ 



C. Doped hole: 
coupling with magnons and orbitons 

We expect that a doped hole does not modify signifi- 
cantly the classical ground state stable for the half-filled 
case (see above). This could play a role in the lightly 
doped regime, but in the case of one hole in the whole 
plane such a modification is negligible and will be ne- 
glected below. Instead, the doped hole may modify its 
neighborhood by its coupling to the excitations of the 
classical ground state — magnons and orbitons — which 
renormalize the hole motion. In order to describe it 
mathematically, we rewrite Ht (see below) and H^s (see 
next section) using similar transformations as performed 
in Sec. HlBl 

First, we rotate spins and pseudospins on sublattice A. 
Next, we express the electron operators in terms of the 
Schwinger bosons introduced in Sec. IIIBI 



1 



P - —fU^h. 



(2.37) 
(2.38) 
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with the same constraints as in Sec. IIIBI Note that the 
factor is not added "ad hoc" but originates from a 
detailed check of the vaUdity of the above equations: It 
should always be present in the case of spin = 1 be- 
cause e.g. when one annihilates two-boson state with the 
/ operator, then a factor appears. Due to this fac- 
tor and the above constraint on the number of bosons 
the high spin projection operators V in Ht arc no longer 
needed (i.e., quantum double exchange^ factors are im- 
plicitly included in this formalism). 

Next, similarly as in Sec. IIIBI we transform the 
Schwinger bosons into the Holstein-Primakoff bosons, 
skip all terms containing more than two bosons, perform 
Fourier transformation for bosons and (additionally) for 
holons here, introduce ak± operators, and finally perform 
Bogoliubov transformation to arrive at the Hamiltonian: 

r 

'^*"2]V ^ P'^:^^(k,ql,q2)^L/^kB(aql+ +aqi-)/3q2A 

k,qi,q2 

-hMy(k,qi,q2)/i[s^kA("qi+ - aqi-)/3q2S 
-hM^(k, qi , q2)^L'^kB(a-qi+ - "tqi- )P^2A 

(2.39) 

where 

M,,(k,qi,q2) = Uqi7fc^,-gi,,-g2^, + «qi7fc^,-g2,, > (2-40) 

with jjL = x,y and k = k— qi — q2; the coefficients 
{Uq^,Vq^} are the standard Bogohubov factors (|2.30p . 

D. Doped hole: the free dispersion 

After performing similar transformations as the ones 
introduced in Sec. Ill CI we obtain that the three-site 
terms, Eq. (|2.10p . lead to the following Hamiltonian for 
the holes 

'^3s = rJ2 {£Bik)hlsh^,B + eAik)hl^h^,^} , (2.41) 

k 

where 

r=ii^J, (2.42) 

and the free dispersion relations for two orbital flavors 
are: 

EAik) = 2cos{2ky), (2.43) 
eB(k) = 2cos(2A::,). (2.44) 

Note, that, we neglected all of the three-site terms which 
lead to the coupling between holes and magnons and/or 
orbitons. This approximation is physically well justified 
since all such terms would be of the order of J/4, i.e., 
much smaller than the terms in Eq. (|2.39p . Sec also 
Appendix [Bl for further discussion. 
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FIG. 2: (Color online) The three-site hopping term r (|2.42p 
(solid line), the effective spin exchange Js (|2.35l) and orbital 
exchange Jo (|2.34|) interaction (dashed and dotted line) for 
increasing Hund's exchange coupling t? (|2.9p . The realistic 
value of = 0.15 for LaVOs (cf. Ref. wW l is indicated by the 
light vertical dotted line. 



E. Polaron model and its parameters 



The above considerations demonstrate that in the 
lightly doped case, when the classical spin and orbital 
ordered ground state present in the half-filled case is ro- 
bust, the strong-coupling model (|2.ip can be reduced to 
an effective model 



n = nt + Hj + n3s, (2.45) 



see Eqs. ((23T1) and ((2?39)) - ((2^ . Actually, this is a 
polaron-typc model with the coupling between fcrmions 
(holes) and bosonic excitations (orbitons and magnons) 
which is relatively easy to solve, see below. The validity 
of the mapping between the two models was thoroughly 
discussed in Ref. d and [H. 

Note, that the original strong-coupling model (|2.ip has 
three parameters {J,r],t}, whereas the effective polaron 
model given by Eq. (j2.45p is more conveniently analyzed 
using instead four parameters {Js, Jo^T^t}, which deter- 
mine the scale of spin and orbital excitations, as well as 
free hole propagation (due to the three-site terms) and 
the vertex function (t), see below. In what follows we 
will use cither one of these two parameter sets (or even 
both of them), depending on the context. Hence, we plot 
in Fig. [2] the functional relation between the parameters 
{J5, Jo,t} on Hund's exchange rj. While the magnon 
energies oc Js decrease with 77, the energy scale of or- 
bitons (X Jo increases rapidly, so the latter excitations 
are expected to play an important role in the realistic 
regime of parameters. 
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III. GREEN'S FUNCTIONS AND SPECTRAL 
FUNCTIONS 

A. The self-consistent Born approximation 

The spectral properties of the hole doped into the 
AF/AO ground state |<i>o) with energy Eq of the strong- 
coupling model Eq. at half-filling follow from the 
Green's functions: 



Gai(k,o.) 



$0 
$0 
$0 
$0 



1 



U} - 




-Eo 




1 




LU - 




-Eo 




1 




Lu ^ 


-H 


-Eo 




1 




LU H 


-H 


-Eo 



flki 


$0^ 


, (3.1) 


flkT 




, (3.2) 


bkl 




, (3.3) 


fokT 


$0^ 


• (3.4) 



However, due to the mapping of the strong-coupling 
model onto the polaron model performed in Sees. IIIBh 
IIIE| it is now convenient to express the above Green's 
hmctions in terms of the polaron Hamiltonian Ti.. This 
requires that we first write down the electron operators 
in terms of the operators used in Eq. p.45p : 



1 



Oki = 



OkT 



'ki 



&kT = 



TV 



N 



N 



(3.5) 
6) 
(3.7) 



TV 



(3.8) 



Second, the ground state |<I>o) is now a physical vacuum 
|0), with respect to the Bogoliubov operators liki and 
the operators /3k, with energy E. Then, we arrive at the 
following relations: 



Ga4(k,C.) = 1(0 



ikA 



u + H-e'"'^^ 



GAA(k, w), 



G,T(k,c.) = 1(0 



tkS 



Lu + n-E^"^^ 



GBB(k,w), 



whereas 



GaT(k,t^)<Gai(k,w), 
Gbi(k,w) < G6|(k,w), 



(3.9) 



(3.10) 



(3.11) 
(3.12) 



M, 



ORBITON 
" MAGNON 




FIG. 3: (Color online) Diagrammatic representation of the 
SCBA equations: top — the Dyson's equation for the 
GyiA(k, tj) and Gflfl(k,cj) Green's functions; bottom — equa- 
tions for the respective self-energies. Note the appearance of 
the two wiggly lines in the calculation of the self-energies com- 
ing from simultaneous orbiton and the magnon excitations. 



and we skip the two latter Green's functions in what 
follows. Note that the above set of equations follows 
from the fact that /3k |0) = and the inequalities are due 
to 



^-0.1, 



(3.13) 



where tiq is the average number of spin deviations in the 
2D ground state |0). Note that below we will eliminate 
the ground state energy E to simplify equations. 

As seen in Section III CI the vertices in the spin-orbital 
model are slightly more complex than for the standard 
spin t-J modeli^ (i) one always finds two boson and two 
holon lines at each vertex (instead of just one boson and 
two holon lines), (ii) the two sublattice structure plays an 
essential role (this resembles the orbital case), and (iii) 
one has two kinds of magnons (which we had to intro- 
duce in order to keep track of the lattice index). Also, 
in the LSW and LOW order all the terms cx Q;;/3j do not 
contribute (the self-energies for them would require alto- 
gether four boson lines instead of just two for cx a-^Pi)^ 
Hence, we largely follow the route proposed for the or- 
bital t2g modeli^ and obtain the following SCBA equa- 
tions for the self-energies (see also Fig. [3|): 



S]^A(k,c.) = 2j^ ^ M2(k,qi,q2) 



4/V 



qi.q2 



X GBB(k - qi - q2,t^ + ./s^qi + Jq) 



2+2 



EBB(k,^) = 2^ ^■^,'(k,qi,q2) 



4/V2 



qi:q2 



X Gaa{^ - qi - q2,'^ + ^st^qi + ^o) , 



(3.14) 



(3.15) 



where the factor 2 in front of each vertex comes from the 
fact that we have two kinds of magnons [-1- and — , see Eq. 
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(I2.28p ] and hence two distinct diagrams. Fortunately, 
this factor cancels with one of the 2 in the denominator 
and we obtain that the coupling constant is {t/\/2)^, i.e., 
we recover this factor 1/^/2 which originates from the 
quantum double exchangei^ The above equations should 
always be supplemented by the Dyson's equations: (Fig. 



where we introduced a factor of 2 in front of the definition 
of the spectral function A-y(k, w) for convenience. 



G_BB(k,t^) - 



1 



- r£A(k) - i;AA(k,a;) 
1 

reB(k) - I]sB(k, w) 



(3.16) 
(3.17) 



B. Analytic structure and quasiparticle pole 



which together form a self-consistent set of equations 
which can be solved numerically by iteration. 

Finally, we can calculate the spectral functions for a 
hole created in a and b orbital: 



^a(k,w) 



- — lim Im Gn i (k, 



i6) 



lim ImGAA(k,tJ + i(5), (3.18) 



Ab(k, w) = lim Im Gfet (k, w 

TT 5^0 
1 



t6) 



- lim ImGi3s(k,cj + i(5), (3.19) 



It occurs that in the case when the three-site terms 
are absent (i.e., at r = 0) one can easily prove two im- 
portant properties of the SCBA equations (|3.14p - (|3.15p : 
(i) the self-energies are k- independent, (ii) the spectral 
functions contain the quasiparticle state for finite value 
of the exchange parameter J (equivalently J5 or Jq). 

First, we show property (i). Since we assumed that 
T = 0, we can rewrite SCBA equations (|3.14p - (|3.15p to- 
gether with Dyson's equations (|3.16p - (|3.17p in the follow- 
ing manner: 



EyiA(k,tj) 
I]sB(k,t^) 



A/2(k,qi,q2) 



2iV2 uj + Jswqi + Jo - SsB(k - qi - q2, w + Js^qi + -'^o) 



2+2 



2N^ ^ uj^ 
qi,q2 



E 



Af2(k, qi,q2) 



Jsi^cix + Jo - S^yi(k - qi - q2, + ^s^qi + Jo) 



(3.20) 
(3.21) 



which after substitution q2 ^ k — q2 in the sums leads to 



E 



/'(qi,q2) 



2iV^ qf^2 "^^^qi + Jo - S_BB(q2 - qi,uj- Js^qi - Jo)' 



24-2 



.9^(qi,q2) 



q7^2 ^ ~^ -^s^qi + Jo - ^aa{(12 - qi, UJ + Js^^q^ + Jo) ' 



(3.22) 
(3.23) 



where we have used a short-hand notation: 

/(qi,q2) = Uqilq2.-qi. + Vcillq2. ^ (3-24) 
5(qi,q2) = Uq^lq2y-qiy + Vciilq2y ■ (3-25) 

Since the self-energies on the right hand side of Eqs. 
(|3.22p ~ (|3.23|) do not depend on k we are allowed to drop 
the momentum dependence of the self-energies. Hence, 
the spectral functions are also momentum independent 
in the case of r = and we recognize that the depen- 
dence of the spectral functions on k may originate only 
from the three-site terms. 

Second, using the dominant pole approximation^ we 



show that the quasiparticle state exists [property (ii)] if 
J is finite (i.e. Js or Jo are finite). Hence, following 
Kane et al.^ we assume that the Green's function can be 
separated into the part containing the pole and the part 
responsible for the incoherent processes: 

GAAii^) = -^^ + G'X1{^), (3.26) 
UJ — uja 

Gbb(w) = ^^ + G*^^(w), (3.27) 

Ll) — UJB 
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where 



aA 



a-B 



and the pole positions: 



UJA 
LOB 



1 

duj \u—lob 



(3.28) 
(3.29) 



(3.30) 
(3.31) 



have to be determined self-consistently following Eqs. 

Next, following Ref. [3, it is straightforward to derive 
the upper bound for the residues {aA^aB}'- 



2j-2 



a-B 



2N^ 



aB < <! i + E5'(qi'q2) 



2iV2 



(Js^^qi + Jo) 



(3.32) 

-1 



(3.33) 



If the sums in the above equations are divergent, then 
the upper bounds for the residues are equal zero and 
the Green's functions do not have the quasiparticle pole. 
Hence, we need to check the behavior for small val- 
ues of the momenta qi. Then cx |qi|, but e.g. 

("qi7<j2„-9l„ +Wqi792j^ ^ HMl q2y ' ^1 ' ^ 1 qiyf ■ ThuS, 

if at least either J5 or Jq is finite, then there are no di- 
vergences in Eqs. (|3.32p - (|3.33|) . Consequently, under the 
same conditions the quasiparticle state exists. 



IV. NUMERICAL RESULTS 

A. Spectral functions 

We calculated the spectral functions yla(k, lj) and 
Ab(k,w) by solving SCBA equations ((3?T4)) - ([3?T5l) on a 
mesh of 16 x 16 k-points. The results for two different val- 
ues of the superexchange constant J = 0.2i and J = Q.Qt 
are shown in Figs. S] and O respectively. Also we as- 
sume the Hund's coupling to be close to a realistic value 
in LaVOa, i.e., -q = 0.15.— We emphasize that such a 
value of the Hund's coupling is not only realistic but to- 
gether with the observed value of J = 0.2t gives a reason- 
able value of the spin-only exchange constant J5 — 0.03t 
(which is in agreement with the observed value of the 
Neel temperature in LaVOg Tn ^ 143 K)^L^ cf. Sec. 
VIB of Ref. m for a more detailed discussion on this 
issue. Thus, the spectral functions shown in Fig. |4] are 
calculated for the realistic values of parameters in LaVOs 
whereas those shown in Fig. [S]will serve as a comparison 
with the regime of large J. 



3 

< 



3 




FIG. 4: The spectral functions derived within the SCBA for 
the hole doped at half-filling to the spin-orbital model (|2.1|) 
into: (a) a orbital, and (b) b orbital. Parameters: J = 0.2t 
and = 0.15 (i.e., Js = 0.03t, Jo = O.OSt, and r = OMt). 
Broadening S = O.Olt and cluster size 16 x 16. 



Let us first discuss the results presented in Fig. [4] The 
quasiparticle peak in the low energy part of the spectrum 
is clearly visible and confirms the analytic predictions 
presented in Sec. IIIIBI However, the quasiparticle state 
has a rather weak ID dispersion: along the direction 
for holes doped into the b orbitals and along the ky di- 
rection for holes doped into the a orbital. Since such 
a dispersion was not present when the three-site terms 
were neglected (see Sec. IIII B[) we can immediately as- 
cribe the onset of the ID dispersion in the spectra to the 
three-site terms. Indeed, these terms are responsible for 
the coherent hole propagation in the ground state with 
AO order and the hole dispersion is dressed by incoherent 
processes. This phenomenon is quite well understood in 
the framework of the orbital-only strong-coupling model 
— we refer to Ref. [l^ and Ref. [l^ for the more detailed 
discussion. 

Besides, the excited states form a ladder-like spectrum, 
also with a small ID dispersion. Similarly to the quasi- 
particle state the excited states resemble qualitatively the 
spectral functions calculated for the purely orbital model, 
see Ref. [l2| and Ref. [ij. However, here we see that the 
spin-orbital spectra are quantitatively different than the 
orbital ones. In particular, they are quite similar to those 
obtained for J = OAt in the purely orbital model (see Fig. 
9 in Ref. [Tsh although this is just a coincidence since all 
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co/t 

FIG. 5: The spectral functions derived within the SCBA for 
the hole doped at half-filling to the spin-orbital model (|2.1|) 
into: (a) a orbital, and (b) b orbital. Parameters; J = 0.6t 
and ri = 0.15 (i.e., Js = 0.09t, Jo = 0.23t, and r = 0.19t). 
Broadening S — O.Olt and cluster size 16 x 16. 



of the exchange constants used in the spin-orbital model 
(J, Js or Jo) are much smaller than 0.4t. 

Looking at the spectral functions shown in Fig. [5] 
we note that the ladder-like spectrum almost disappears 
and the quasiparticle dispersion increases to ^ OAt when 
the value of the superexchange constant is increased to 
J = 0.6t. However, the qualitative features in the dis- 
persion of the quasiparticle peak stay mostly unchanged 
and the ID character of the dispersion relation is pre- 
served. Quantitatively, both the quasiparticle spectral 
weight and the bandwidth increase quite drastically. 

To conclude, the spin-orbital spectral functions form 
ladder-like spectra with a small ID dispersion and have 
many similarities with the purely orbital spectra of the 
t2g modelJ^ Still, however, there are few relatively im- 
portant differences with the orbital model which suggest 
that the spin-orbital case is more complex and the quasi- 
particle behavior more subtle than for the purely orbital 
model. For the understanding of this problem we refer to 
Sec. |V]whercas in the next section wc further investigate 
the properties of the quasiparticle arising in the spectral 
functions of the spin-orbital model. 



B. Quasiparticle properties 

In this section we analyze in more detail the quasiparti- 
cle properties of the spectral functions of the spin-orbital 
model, cf. Fig. [S) In addition, we compare these prop- 
erties with those calculated for the spin strong-coupling 
model of Ref. [13 and the orbital strong-coupling model 
of Ref. [H. 

By looking at the results for the spin-orbital model 
in Fig. [S] we see that: (i) the quasiparticle bandwidth 
W is strongly renormalized from its respective free value 
(which is = 8r as it may only originate from the three- 
site term dispersion relation, sec Sees. IIII BlllV^ and is 
proportional to for small J {J 0.6t) and to J in 
the regime of large J {J ^ 0.6t); (ii) the quasiparticle 
spectral weight aqp grows considerably with increasing 
J; (iii) the pseudogap A (the energy distance between 
the quasiparticle state and the first excited state) exists 
and roughly scales as t{J/t)'^^^ although for larger values 
of superexchange J some deviations from this law were 
observed. 

Thus, on the one hand, we see that the quasiparti- 
cle properties of the spin-orbital model are qualitatively 
distinct than those of the spin model. Wc should stress 
that although the spin model used here is the 5=1/2 
t-J model with the three-site terms-- but we expect sim- 
ilar generic behavior in the case of spin 5 = 1 (see also 
Sec. IV Ap . Hence, as there is no qualitative agreement 
with the model for spin 5=1/2 there would neither be 
a qualitative agreement with the model for spin 5=1. 

On the other hand, we note remarkable similarities be- 
tween the present model and the orbital model (see also 
Fig. 0; the generic behavior of the bandwidth, the spec- 
tral weight, and of the pseudogap seems to be almost 
identical in both models. There are, however, two im- 
portant differences. First, the t{J/t)'^/^ law docs not de- 
scribe the behavior of the spin-orbital pseudogap so well 
as in the orbital case. Second, if wc assume that the 
spin-orbital model is just qualitatively similar to the or- 
bital model, then it should be possible to rescale all the 
quasiparticle properties with such an effective value of 
the superexchange J that they coincide with the results 
for the orbital model. This, however, is not possible. For 
example, from the fits to the t{J/t)'^/^ law we can de- 
duce that such an effective value of the superexchange 
would be Jeff = (aso/ao)^^^ ~ 0.79 J, where aso = 1-94 
{ao = 1.66) is the fitted coefficient which multiplies the 
t{J/t)'^/^ law in the spin-orbital (orbital) case. This sin- 
gle parameter does not suffice, however, as a similar fit 
to the quasiparticle spectral weight would require such a 
Jeff to be bigger than J. 

In conclusion, the properties of the quasiparticle state 
in the spin-orbital model resemble those found for the 
quasiparticle in the purely orbital t2g modelJ^ii^ How- 
ever, a more detailed analysis presented in this section 
reveals that such a correspondence is rather " superficial" 
and that few of the qualitative features of both models 
are different. We refer to the next section for the thor- 
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J/t 

FIG. 6: (Color online) Quasiparticle properties for increasing 
J for the spin-orbital model discussed in the present paper 
with the Hund's coupling rj — 0.15 (solid line), compared with 
the spin model of Ref. [sj (dashed line), and with the orbital 
model of Ref. [l^ (dotted line). Different panels show: (a) 
the quasiparticle bandwidth W, (b) the quasiparticle spectral 
weight aqp averaged over the 2D Brillouin zone, and (c) the 
pseudogap A (energy distance between the quasiparticle state 
and the first excited state) at k = {ir /2,n/2). Note that A is 
not shown for the spin model as it cannot be defined there. 
The two light solid lines on panel (c) show the (J/t)^^^ curves 
fitted to the data at low values of J for the orbital and the 
spin-orbital model. 



ough understanding of this phenomenon. 



V. ORIGIN OF SPIN-ORBITAL POLARONS 

A. Comparison with spin and orbital polarons 

The main task of this and the following two sections 
is to understand the spectral functions and the quasipar- 
ticle properties of the spin-orbital model. In particular, 
we not only want to understand the peculiar similari- 
ties or rather small differences between the spin-orbital 
model and its purely orbital counterpart, which were dis- 
cussed in the last two sections (see Sees. IIV Al and IIVB|) . 
but also we want to understand why the spin physics 



seems to be "hidden" in the present spin-orbital system. 
Note that for the sake of simplicity in this section we will 
entirely neglect the three-site terms (|2.10p in the spin- 
orbital model (|2.ip (or in other words we will put r = 0). 
This is motivated by the fact that the role of the three- 
site terms is solely to provide dispersion to the spectra 
(see Sees. IIII BI IIV Bp and the origin and character of 
these dispersive features is quite well understood by now, 
see. Ref. M 

We start the analysis of spectral functions by intro- 
ducing two artificial toy models in order to extract the 
features related to spin and orbital excitations separately, 
whose results will be later compared to those of the t-J 
spin-orbital model. First, we define the following t-Js 
toy-spin model 

Hs=-tJ2 P(4cj<x + H.c.)P + J5 5^ Si- Sj, (5.1) 

(ij) 

where spin operators {Si} stand for S ~ I spins, the 
constrained operators, C;^ = C;g.(l —ni^), exclude double 
occupancies from the Hilbert space, and the operators 
V project onto the high spin states. Note that here the 
superexchange energy scale is Js [see Eq. (|2.35p ] and not 
J. Hence, it is defined in such a way that it mimics the 
influence of AO order on the spin subsystem. On the 
other hand, the kinetic energy is blind to the AO here 
and Eq. (|2.2p reduces to the kinetic part of Eq. (|5.ip 
only if the orbitals form an orbital liquid state. This is 
an obvious logical inconsistency but our aim here is to 
see what happens when the joint spin-orbital dynamics 
in the kinetic energy is entirely neglected. It is also the 
reason why we call model (jS.ip the toy-spin model. 
Second, we define the following t-Jo toy-orbital model 

i (ij) 

(5.2) 

where pseudospin T = 1/2, and the constrained opera- 
tors 6; = 6; (1 — nia) and a; = a; (1 — nib) exclude double 
occupancies. Similarly as in the toy-spin model defined 
above, the superexchange energy scale is not J but Jo 
[see Eq. (|2.34p ] which mimics the influence of spin system 
with the AF order on the orbital one. Also, the kinetic 
energy is blind now to the AF order and Eq. (|2.2p reduces 
to the kinetic part of Eq. (|5.2[1 only if the spins form the 
FM order. This again is logically inconsistent, but we 
wish to consider this situation for the same reasons as 
for the spins, see above. 

Next, we solve these models using the SCBA method 
in the case of the one hole doped into the AF (AO) state 
for the toy-spin (toy-orbital) model. We do not show here 
the respective SCBA equations as these are straightfor- 
ward and follow from those written in Ref. i (Ref. IH) m 
the case of the spin (orbital) model. One only has to re- 
member to substitute in the respective SCBA equations 
J — > — J5, S* — !• 1 and (due to the quantum double ex- 
change factor) also t t/y/2 in Ref. in the spin case. 
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FIG. 7: (Color online) Spectral function along the F — M 
direction of the Brillouin zone for: the spin-orbital model (top 
panels), toy-spin model (middle panels), toy-orbital model 
(bottom panels). Parameters: J = 0.2t, r] = Q (i.e., Js = 
O.OSf, Jo = 0) on the left panels whereas J = 0.2t, 77 = 0.15 
(i.e., Js = 0.03t, Jo = O.OSt) on the right panels; in both 
cases T = 0, see text. Note that r = implies^ Aa{^,Lo) — 
^6(k,w) = A(k,uj). Broadening 5 = O.Olt [5 = 0.02i on panel 
(f)] and cluster size 16 x 16. 



FIG. 8: (Color online) Spectral function along the T — M 
direction of the Brillouin zone for: the spin-orbital model (top 
panels), toy-spin model (middle panels), toy-orbital model 
(bottom panels). Parameters: J = 0.6t, rj — (i.e., Js = 
0.15t, Jo = 0,) on the left panels, whereas J — 0.6t, 77 = 0.15 
(i.e., Js = 0.09t, Jo = 0.23t) on the right panels; in both 
cases r = 0, see text. Note that r = implies^ Aa(k,aj) = 
Ab(k,uj) = A{'k,u}). Broadening 5 = O.Olt [5 = 0.02t on panel 
(f)] and cluster size 16 x 16. 



and J ^ — Jo, Eo ^ 0, and r ^ in Rcf. [l^ in the 
orbital case. We then calculate the respective spectral 
functions on a mesh of 16 x 16 k-points. Note that since 
T = the spectral functions for both orbital flavors are 
equal, i.e. Aa{li,uj) = Ab{ls.,uj) = A(k, a;). 

Finally, we compare the results obtained for the above 
toy models with those obtained for the spin-orbital t-J 
model, Eq. (|2.2p and (|2.5p . We show the results for two 
different values of J = 0.2t and 0.6i, see Figs. [7]and[51 In 
addition, we calculate the results for two different values 
of the Hund's coupling ij = and rj ~ 0.15 (see left and 
right panels of each figure). 

Let us first look at the physical regime of finite 77 = 0.15 
and J = 0.2t, see Figs. [7i;d)-17i;f). On the one hand, the 
spin-orbital spectral function [panel (d)] resembles qual- 
itatively the ladder spectrum found in the orbital model 
[panel (f)] although the quantitative comparison reveals 
strong differences between the two models. On the other 
hand, the spin-orbital spectral function is entirely differ- 
ent from the k-depcndent spin spectral function [panel 
(e)]. One also finds a somewhat similar behavior for the 
case of 77 = 0.15 and J = OM, see Figs. [5i;d)-[Hi;f). Here, 



however, the spin-orbital spectrum is much more differ- 
ent than the orbital spectrum, so the rcnormalization by 
the spin part is stronger. 

Even more inquiring behavior is found in the nonphys- 
ical regime of 7; = (which, however, is an interesting 
mathematical limit) 1^ Then neither of the panels shown 
in Figs. [7]^a)~[7Jc) or Figs. [5{a)-[5]Jc) is similar to each 
other. This means that even the orbital model is entirely 
different in this regime than the spin-orbital model. This 
is because in this limit the hole moves in the orbital model 
incoherently as Jo = for 77 = 0, see e.g. Fig. Ulc). How- 
ever, apparently in the spin-orbital model with Jo = 
and small but finite J5 the hole moves in string-like po- 
tential, see e.g. Fig. Et^a). This means that the onset of 
the ladder-like spectrum in the spin-orbital model in this 
regime cannot be explained easily in terms of the purely 
orbital model. 

In addition, one sees that whereas for different values 
of 77 but the same values of J one gets rather similar 
spin-orbital spectra, the spectra found for the toy-orbital 
model are different. On the contrary, for different values 
of J and the same values of 77 the spin-orbital spectra are 
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rather distinct but the spin spectra do not change much. 
This is another argument which suggests that neither 
the toy-spin model nor the toy-orbital model can explain 
alone the features found in the spin-orbital spectra. 

To conclude, we note that the joint spin-orbital dy- 
namics in the kinetic part of the spin-orbital model plays 
a significant role for the coherent hole motion. When 
spin and orbital degrees of freedom are separated from 
each other (disentangled) , the purely spin or orbital mod- 
els cannot reproduce the spectral function found for the 
spin-orbital t-J model. Hence, indeed the spin-orbital 
spectral functions resemble the orbital ones only super- 
ficially and it is the peculiar interplay of the spins and 
orbitals, studied in the next section, which leads to the 
calculated spectra. 



B. Spin-only string excitations at = 

In this section we attempt^ to understand the spin- 
orbital spectra in the limit of = by assuming that 
the spins 5 = 1 are purely classical and interact by 
Ising superexchange. Hence, we skip the spin-fiip terms 
(X S'^S7 in Hamiltonian (|2.5p . rewrite SCBA equations 
(|3.14p - p.l5p in this case, and finally try to compare the 
spectral functions calculated in this regime with the ones 
obtained for the full model, given by Eq. (|2.1|) . In addi- 
tion, we not only assume 77 = (which implies Jo = in 
particular) but also t = (similarly as in the previous 
section). In the next section we discuss the impact of the 
finite value of these parameters on the results obtained 
under these assumptions. 

Since Wq^ = 1, Wq^ = 0, and Wq^ = zS for the 
Ising superexchange^^ we can rewrite self-energy equa- 
tions (IS2ni)-(E2Il): 



i;AA(k,w) 



zH'' 



E 



7, 



27V2 ^ UJ + JszS -Y.BB{cii,^ + JszS)' 

(5.3) 

^Si3(k,t^) = 2iV^ E^ w + Jsz5 - Y^aMi ,lu + JszS)' 

(5.4) 

where we already substituted qi — > k — qi — q2 in the 
sums. Then the above self-consistent equations are mo- 
mentum independent and we obtain: 



Sbb('^) 



-) - 

^/2/ uj- 



z/2 



JszS - ^Bsi^ 
z/2 



JszS) 



JszS - Yaa{uj + JszS) 



(5.5) 



(5.6) 



since 1/A^ Y^qi ^'qi^ ~ V(2-2)j where 1/ = x,y. 

We have solved Eqs. (|5.5p - (|5.6p self-consistently. The 
respective spectral functions are shown in Fig. [5] for 
J = 0.2t (i.e., Js = OMt) and J = 0.6t (i.e., Js = 0.15t). 




co/t 

FIG. 9: (Color online) Comparison between the spectral func- 
tions of the classical limit (dashed lines) and full quantum 
version (solid lines) of the spin-orbital model at a single point 
in the Brillouin zone (the spectra are k-independent) for: (a) 
J = 0.2t and 77 = (i.e., Js = O.OSt, Jo = 0), and (b) J = 0.6t 
and 77 = (i.e., Js = 0.15t, Jo = 0); in both cases r = 0, 
see text. Note that r = impliesi^ Aa{io) = Ai,{ll)) = A{lu). 
Broadening S = O.Olt {S = O.OlSt) in panel (a) [(b)] and clus- 
ter size 16 X 16. 



As expected, one obtains a typical ladder spectrum in 
the entire regime of J. However, we see that the results 
resemble those obtained for the full spin-orbital model 
(|2.ip with = and r = 0: Although the spectrum of the 
full spin-orbital model contains some incoherent part, the 
ladder peaks of the full spin-orbital model and of its clas- 
sical version (with Ising spin superexchange) almost coin- 
cide. In addition, the incoherent bandwidth in the J ~ 
limit is W = At in the classical case whereas it is only 
shghtly reduced in the quantum model {W ~ 3.7t)^ For 
finite J this results in the small shift of the peaks in the 
full spin-orbital model with respect to its classical coun- 
terpart. Altogether, this suggests that the classical and 
the full (quantum) versions of the spin-orbital models are 
to a large extent equivalent. 

We note that Eqs. (|5.5p - (|5.6p are almost identical to 
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the SCBA equations for the hole moving in the S = 1/2 
spin Ising model [cf. Eq. (20) in Ref. The only dif- 
ferences are: the self-consistent dependence of the self- 
energies on two different sublattices, the reduction of the 
nearest neighbors by a factor 1/2 (in the numerator), 
the reduction factor 1/V2 for the hopping, and the in- 
crease of the magnon excitation energy by a factor of 2. 
Whereas the first two imply the zig-zag hole motion in 
the ordered state^^ the two others merely mean that the 
hole moves in the spin S = \ system, cf. similar factors 
for the triplet hole in the t-J modeli^ This demonstrates 
that the hole motion in the full spin-orbital model with 

= and r = can be quite well approximated by the 
zig-zag hole motion in the S = \ spin Ising model. 

The whole analysis presented above leads us to the 
conclusion that in the limit of Jo = and t = the 
hole moves in the spin and orbitally ordered plane in the 
following way: (i) the orbitals force the hole to move 
along the zig-zag paths even in the limit of Jq = 0, (ii) 
the orbitals force the hole to retrace its path again even 
in the limit of Jq — > — this resembles the situation 
discussed by Brinkman and Rice^ for the spin system, 
where the hole in the Ising model with J = always 
has to retrace its path, (iii) the coherent hole motion by 
the coupling to the spin fluctuations cannot occur in the 
ground state as then the hole would not retrace its path, 
so (iv) instead the hole creates strings in the spin sector 
which are erased when the hole retraces its path. Thus, 
we note here a complex interplay of spins and orbitals. 
In particular, due to point (ii) the orbitals force the spins 
to effectively act on the hole as the classical objects. 



C. Composite string excitations at > 

One may wonder how to extend the above understand- 
ing of the role of spin and orbital excitations in the forma- 
tion of spin-orbital polarons to the case of finite values 
of orbital exchange interaction Jo > at 77 > (see 
Fig. [2|), or in the presence of three-site terms at finite 
T. On the one hand, including the nonzero value of Jo 
merely leads to the substitution of JszS — s- JszS + Jo 
in Eqs. (|5.3p - (|5.4p and consequently also in Eqs. (|5.5p - 
(|5.6p ; this increases the string- like potential acting on the 
hole, coming in this case also from the orbital sector. On 
the other hand, including the three-site term r > re- 
sults in the energy shift uj — *■ cj-|-e^(k) and uj uj+eBOi-) 
in Eqs. (|5.3p - (|5.4p which results in the k-dependence of 
the spectral functions, see Eqs. ()5.5p - (|5.6p . However, 
we can still solve the model using the SCBA. The re- 
sults (not shown) resemble those found in Figs. |4] and 
[SJ it is again only the incoherent part which is slightly 
enhanced in the full spin-orbital model (j2.ip whereas in 
its classical counterpart it is suppressed. Furthermore, 
all of the quasiparticle properties of the full spin-orbital 
model shown in Fig. [5] are almost perfectly reproduced 
by the classical spin-orbital model (not shown) — with 
the only slight discrepancy occurring in the region of the 



slight deviation from the t{J/t)'^/^ law in the pseudogap 
of the full spin-orbital model. 

Finally, there remains just one subtle issue: how to 
explain the appearance of the small momentum indepen- 
dent incoherent part in the high-energy part of the spec- 
trum for the quantum spin-orbital model, cf. Fig. [5] in 
the case oi rj = [a similar incoherent dome is visible 
in the spectrum for finite 77 (not shown)]. This can be 
understood in the following way: although the hole has 
to return to the original site (due to the orbitals) the 
magnons present in the excited states can travel freely 
in the system. Hence, the energies of the excited states 
can no longer be classified merely by the length of the 
retraceable paths (as it would be the case in the classi- 
cal model with no dispersion in the magnon spectrum). 
This results in the small incoherent spectrum which sur- 
rounds each of the peak of the ladder spectrum visible in 
for example Fig. [5) Furthermore, this incoherent spec- 
trum grows with increasing J5 as then the velocity of the 
magnons increases. 

To conclude, we note that the quantum spin fiuctua- 
tions are to a large extent suppressed in the spin-orbital 
model by the simultaneous coupling of the hole to both 
spin and orbital excitations. In particular, they do not 
affect the quasiparticle state and merely add as a small 
incoherent spectrum in the incoherent high-energy part 
of the ladder spectrum. This is due to the classical char- 
acter of the orbitals which confine the hole motion and 
prohibit its coherent motion by the coupling to the quan- 
tum spin fluctuations. On the other hand, the hole still 
couples to the spin degrees of freedom, mostly, in a classi- 
cal way, i.e., by generating string potential due to defects 
created by hole motion. Thus, the string which acts on 
the hole moving in the plane with AO and AF order is of 
the composite orbital and spin character. This not only 
explains the peculiar correspondence between the orbital 
and spin-orbital model but also explains that the spins 
" do not hide behind the orbitals" but play an active role 
in the lightly doped spin-orbital system. 



VI. CONCLUSIONS 

Let us now summarize the results of the present study. 
We studied the motion of single hole doped into plane of 
the Mott insulator LaVOa, with coexisting AF and AO 
order^ shown in Fig. 1, which violates the Goodenough- 
Kanamori rules. The framework to describe the hole 
motion in this situation is the respective t-J model 
with spin-orbital supercxchangCr^iSl supplemented by 
the three-site terms required here^ (which were derived 
in the Appendices) and we used this model as the start- 
ing point of our analysis. Next, we reduced this model 
to the polaron-type model and showed that the hole cou- 
ples in this case simultaneously to the collective excita- 
tions of both the AF state (magnons) and the AO state 
(orbitons). The problem of flnding the spectral function 
was formulated and solved within the SCBA method. It 



15 



should be emphasized that, to the best of our knowl- 
edge, the problem of a hole coupled simultaneously to the 
magnons and orbitons has never been addressed before.— 
We would like to note that the presented spectral func- 
tions, such as the ones shown in Fig. 21 predict the 
main features to be measured in future photoemission 
experiments on the cleaved LaVOa sample. The reader 
may wonder whether (apart from the matrix-elements 
effects responsible for certain redistribution of spectral 
intensity) any other processes, such as for example the 
elcctron-phonon interaction, would affect hole motion to 
such an extent that the spectral functions calculated here 
would change qualitatively. Although we have not made 
any calculations for such a more complex case so far, we 
suggest that they would only enhance the ladder spec- 
trum obtained here, since typically the studied mecha- 
nisms are only responsible for further localization of the 
hole.i^ 

Next, it was calculated both analytically using the 
dominant pole approximation and numerically using the 
SCBA equations that the spectral functions contain a 
stable quasiparticle peak, provided the value of the su- 
perexchange J was finite. We emphasize that this result 
was non-trivial as it was not clear a priori whether the 
coupling between a hole and two excitations would lead 
to a stable quasiparticle state — e.g. the coupling be- 
tween the hole and two magnons does not lead to the 
stable quasiparticle peak, cf. Ref. [H and Appendix |B] 
However, since the orbitons are massive excitations, the 
hole does not scatter strongly on them and the quasipar- 
ticle solution survives. 

Furthermore, we studied in detail the properties of the 
quasiparticle states. In particular, we looked at the differ- 
ences between the well-known spin^ or orbita &^^i^^ po- 
larons and the obtained here spin-orbital polarons. We 
checked that all of the typical quasiparticle properties 
of the spin-orbital polarons such as the bandwidth, the 
quasiparticle spectral weight, and the pseudogap (the dis- 
tance between the quasiparticle peak and the next ex- 
cited state) are qualitatively similar to those of the t2g 
orbital polarons studied recently i^^i^^ and it is the string 
picture which dominates in the quasiparticle properties. 
For example, the bandwidth scales as oc as a result 
of the renormalization of the three-site terms oc J, simi- 
larly as in the purely orbital t2g modelJ^ii^ We suggest 
that the occurrence of the small but still finite bandwidth 
confirms the idea of the absence of hole confinement in 
transition metal oxides with orbital degeneracy presented 
in Ref. [H. 

Finally, a detailed investigation lead us to the conclu- 
sion that the spin-orbital polarons arc microscopically 
much more complex than either spin or orbital polarons, 
and resemble the orbital polarons only superficially. Ac- 
tually, we have shown that the spin degrees of freedom 
also play a significant role in the formation of the spin- 
orbital polarons although their dynamics is quenched and 
they arc forced by the orbitals to act on the hole as the 
classical Ising spins. This happens because the orbitals 



confine the hole motion by forcing the hole to retrace its 
path which means that the hole motion by its coupling 
to the quantum spin fluctuations is prohibited. Thus, the 
string-like potential which acts on the hole is induced by 
the orbitals although it has a joint spin-orbital character. 
Furthermore, it occurred that it is only in the excited 
part of the spectrum that the quantum spin fluctuations 
contribute and are responsible for a very small incoherent 
dome in the spectral function. 

There is also another important feature of this orbitally 
induced string formation: it could be understood as a 
topological effect. This is due to the fact that it happens 
even if the energy of the orbital excitations is turned to 
zero, i.e., when the hole moves in the orbital sector inco- 
herently. Hence, the mere presence of orbitals is sufficient 
to obtain the (almost) classical ladder spectrum of a hole 
doped into the AF and AO ordered state. This result, 
in connection with the fact that the mother-compound 
of the superconducting iron-pnictides shows a variety of 
spin-orbital phenomena,— suggests that further investi- 
gation of the hole propagation in spin-orbital systems is 
a fascinating subject for future studies. 

Having explained all the main goals of the paper, posed 
in the sixth paragraph of Sec. HI we can now try to ex- 
plain the paradox mentioned in the second part of the 
Introduction. Actually, in the last two paragraphs above 
we noted that the spin-orbital polarons resemble superfi- 
cially the orbital polarons but their nature is of joint spin 
and orbital character. On the condition that this result 
could be also extended to the case of the finite doping we 
immediately solve the paradox that experimentally the 
orbitals play a significant role in the doped cubic vana- 
dates, although theoretically we should take into account 
both spin and orbital degrees of freedom. Certainly, any 
straightforward extension of the present one-hole result 
to the regime of finite doping has to be rather specu- 
lative, so further theoretical studies on the doped cubic 
vanadates are needed to verify our explanation of this 
paradox. 
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APPENDIX A: FREE HOLE MOTION DUE TO 
THE THREE-SITE TERMS 

In what follows we derive the three-site terms which 
would contribute to such a motion of a single hole that 
does not disturb the AO and AF order in the ground 
state (as the one described in Sec. IIIB|) . We start our 
analysis by looking at the possibility of the free hole mo- 
tion along the a direction, see Fig. [TOl We would like 
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a orbital 
b orbital 

c orbital 



FIG. 10: (Color online) Three neighboring sites along the a 
direction with such an AO and AF order that the electron on 
site i+a can move over the intermediate site i via the superex- 
change process oc t — J/4 to site i — a without disturbing the 
AO and AF order. 



in Eq. (jAip . then we are allowed to write niai = 1 and 
drop out this operator. Obviously, if we assumed a =1 
or that i Cz B, then we would not get any contribution. 
Thus, for some particular choices of cr and the sublattice 
indices Eq. (|A1[) would describe the free hole motion, 
whereas in some other cases the hopping is blocked and 
it would not contribute at all. While writing down Eq. 
(|2.4ip in Sec. HID I we took care of this problem. 



to move the electron in the b orbital (the one in the a 
orbital does not hop along this direction) from the right 
site to the left site over the middle site in such a way 
that the spin and orbital order present before the pro- 
cess stay intact. This confines the choice of the possible 
high energy intermediate configurations at the mid- 
dle site to the states: (i) the total spin |3/2,— 1/2) 
state with energy U — 3Jh, and the total |l/2, —1/2) spin 
states with (ii) ^E^9 and (iii) '^E^e electronic configura- 
tion, both with energy U. Note that all the intermediate 
states with orbital singlets are excluded as they would 
require orbital excitations. Thus, we arrive at the follow- 
ing contribution to the free hole motion which arises due 
to the three-site term processes along the a direction 



1 1 



3 1-3?/ 



2\ ~. 
-] —Vb 
3/ U 



(Al) 



A similar consideration but for the processes along the b 
direction yields 
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3 1 - 3r/ 
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+h,a 



(A2) 



However, the processes around "the corner" such as 
for instance, first the hopping of an electron along the b 
direction and then along the a direction in the next step, 
would not contribute to the free motion. This is because 
an electron with a particular (a or b) orbital flavor can 
hop only along a single cubic direction in the ab plane, 
and thus we would have to interchange the hopping elec- 
trons at the intermediate high energy site, which would 
lead to an orbital excitation. Hence, after adding the 
sums over all sites i and spins a and the conjugate terms 
to Eqs. ((Il|) - (|J2|) we indeed end up with Eq. (|2.10p 
presented in Sec. Ill Al 

Finally, one may wonder how Eqs. (jAip - (|A2p could 
contribute to the free hole motion since they contain 
four electron operators. However, the number operators 
which stand in the middle of these equations could be 
dropped out for the AO and AF order assumed here. 
More precisely, let us e.g. assume that the intermediate 
site belongs to the sublattice A and concentrate on Eq. 
(|A1[) . Since we assumed in Sec. Ill Bl that in this sub- 
lattice the electrons have spin down and are located in 
the a orbital thus if, in addition, we assume that a =t 



APPENDIX B: COUPLING WITH TWO 
MAGNONS IN THE THREE-SITE TERMS 



A careful analysis leads to the conclusion that all the 
three-site terms which do not contribute to the free hole 
motion would lead to the coupling between a hole and 
either two magnons or two orbitons. On the one hand, 
since orbitons arc local excitations, see Eq. (|2.33p . the 
latter contribution would only slightly enhance the string 
potential in the present model and consequently the spec- 
tral functions would bear even more signatures of the 
ladder spectrum. However, this effect is expected to be 
quantitatively very small as the three-site terms would 
contribute to the vertex as oc (J/4)^, see below, whereas 
the magnon-orbiton vertices considered in the paper are 
of the order of t^. Thus, we can safely neglect these for- 
mer terms. 

On the other hand, neglecting the terms which would 
lead to the interaction between a hole and two magnons 
is a priori not justified. We could imagine that it might 
lead to the hole motion by coupling to the quantum spin 
fluctuations - similarly as in the standard spin case with 
the coupling between a hole and a single magnoni^ Thus, 
we investigated this problem in detail by: (i) deriving the 
respective three-site terms, and (ii) performing all the 
transformations as in Sec lII CI which lead to the Hamil- 
tonian written in the polaron representation. Since all 
these calculations are relatively tedious, we do not write 
down explicitly all the steps but merely present the final 
Hamiltonian which describes the coupling between a hole 
and two magnons 

k,qi,q2 

^ (ctqi+ '^q2+ '^qi+ '-^q^- '^qi- '^q2- '^qi- *^q2+ ) 
+ V'i^(k, qi,q2)/ij^^/iks(aqi + '^q2+ +Q^qi+Q^q2- 

- aqi_ (5q2_ - "qi- aq2+ ) "f" H.C. } 



iT^2^ ^ {l^2,(k,qi,q2)/iL/^kA 



k,qi,q2 



1 Jz 1 

X (<5qi+<5q2+ + '5qi+<5q2_ + ttqi^ttqa^ -|- aqi_aq2+) 

+ Q^q2+ Q^qi+Q^q2- 
+ aqi_(5q2_ - (5qi_Q!q2^) 4- H.C.} , (Bl) 
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where 



following terms: 



Vi^(k,qi,q2) = -Uq.Uq^ cos{2kf, - 2qif, - q2f,) 



+ 2^11^12 cos(2fc^ - q2^) 



(B2) 



V'2,.(k,qi,q2) = -UqiWqj cos(2fc^ - 2gi^) , (B3) 



with all the symbols defined as those in Eq. (|2.39p . Note 
that we neglected above all the terms of the type oc a'^o; 
as they would lead to the four boson line diagrams (the 
self-energies with two magnon lines are very small, see 
below, and hence the self-energies with four magnon di- 
agrams would be even smaller). 

Next, we implement the above derived processes into 
the SCBA method and obtain the following equations for 
the additional self-energies: 

^aa{K^) = ^TTJ Yl {['?^iy(k,qi,q2)-V'2y(k,qi,q2)]^ 



qi.q2 



X GA^Ck- qi- q2,cj+ Jstjq^+ Jswq^) 

+ [?7V^ij,(k,qi,q2) + l/2a(k,qi,q2)]^ 

X Gaa{^- qi - q2, ^st^qi + ^/swqj} , (B4) 
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2iV2 



^ { (k, qi , q2 ) - ^2. (k, qi , q2 )] ^ 



qi.q2 



X G'sB(k- qi- q2, w-l- Jsti^qi+ Jsi^qa) 
-I- [77yi:,(k,qi,q2) + y22;(k,qi,q2)]^ 

X Gbb{^- qi- q2, W+ JsWqi+ Jst^qa)} ,(B5) 



where 
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1 J 
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This requires that we substitute for the self-energies the 



I]ylA(k, 

EBB(k,t^) 
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{k,Lo) + ^AA{Kio), (B7) 
(k,uj)-f EBB(k,w), (B8) 



in the Dyson's equations (|3.16p - (|3.17p which changes the 
SCBA equations. 

Finally, we solve the modified SCBA equations on a 
mesh of 16 x 16 points. It occurs that the spectral 
functions obtained with the additional self-energies Eqs. 
(IB4[) - (jB5p are virtually similar to those obtained with- 
out them (not shown). The only small difference is a 
negligible enhancement of the incoherent part. This can 
be understood in the following way. First, the added 
contributions scale as (J/4)"' and are indeed very small 
as J < i. Second, the vertices, given by Eqs. (|B2p - 
(|B3p . are singular at e.g. qi = (0,0) and q2 — (0,0) 
points. Therefore, the self-energies associated with these 
vertices could only contribute to the incoherent part of 
the spectrum [the divergent vertices lead to the diver- 
gences in denominators of Eqs. (|3.32p - (|3.33p yielding the 
upper bounds for the quasiparticle residues in the domi- 
nant pole approximation] . The physical interpretation of 
this phenomenon is as follows: (i) the hole produces two 
magnon excitation at each step it moves forward, (ii) the 
hole always moves by two sites in a single step (as the 
three-site terms lead to the next-nearest neighbor hop- 
ping), and (iii) magnons "travel" in the system and cure 
the defects created by the hole at a "velocity" one site 
per each step. Thus, one might argue that the magnons 
are not "fast" enough to erase the defects created by the 
hole. 

In conclusion, neglecting the three-site terms which do 
not lead to the free hole motion is entirely justified. Such 
processes are not only quantitatively small but also they 
do not change the physics qualitatively. 
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